Equation of state for Entanglement in a Fermi gas 
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Entanglement distance is the maximal separation between two entangled electrons in a degenerate 
electron gas. Beyond that distance, all entanglement disappears. We relate entanglement distance 
to degeneracy pressure both for extreme relativistic and non-relativistic systems, and estimate the 
entanglement distance in a white dwarf. Treating entanglement as a thermodynamical quantity, 
we relate the entropy of formation and concurrence to relative electron distance, pressure, and 
temperature, to form a new equation of state for entanglement. 
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The relation between thermodynamics and entangle- 
ment is of fundamental importance, and has been studied 
by many authors in the recent years 0, 0> & IS El ■ The 
second law of thermodynamics states that entropy S can 
be regarded as a measure of disorder. It is a macroscopic 
variable and can be related to other state variables like 
pressure and temperature. In quantum mechanics, the 
entropy of formation is a measure of entanglement 0. 
Entanglement is a purely quantum mechanical feature 
and the entropy of formation quantifies how much entan- 
glement is present in a physical system. Entanglement, 
which is usually associated with the microscopic world, 
has recently been shown to also become relevant on 
macroscopic scales 0. Brukner and Vedral showed 
that thermodynamical quantities can also serve as 
entanglement witnesses. All of this suggests that one 
may treat entanglement itself as a thermodynamical 
quantity, and relate it to other macroscopic variables. 
This is the principal motivation of the present work. 

The main system studied in this letter is a Fermi 
gas, i.e. a system of noninteracting fermions. Our treat- 
ment can be extended to other systems, e.g. interacting 
electrons in a superconductor. Entanglement at zero 
temperature for noninteracting fermions has already 
been studied 0, [^j . It was found that all entanglement 
vanishes if the relative distance r between electrons is 
greater then the entanglement distance r e oc 1/kp, where 
kp is the Fermi momentum. In this case, entanglement 
is purely due to particle statistics and not due to any 
physical interaction. At zero temperature, the electron 
gas is in its lowest energy configuration. All energy 
states are occupied up to the Fermi energy ep. The only 
pressure the system exerts is the degeneracy pressure 
P(n) (due to the Pauli exclusion principle), which 
solely depends on its density n. We now want to relate 
this pressure to entanglement. For the non-relativistic 
case, the expression for pressure remains valid if the 
temperature is low compared to the Fermi temperature 



Tp = ep/ks, where ks is the Boltzmann constant. This 
condition is met if the density is sufficiently high, at 
which the fermion gas behaves as an ideal gas. Since the 
particles in an ideal gas are not interacting, we can use 
the same treatment as in 0, llflj . provided the density 
of the fermion gas is high enough. For the extreme 
relativistic case, the expression for pressure is valid for 
all temperatures. 

We will first relate r e to the degeneracy pressure 
P. We then estimate r e in white dwarfs. Third, we 
express entropy of formation and concurrence, which are 
measures of entanglement, as functions of degeneracy 
pressure, relative distance between the electrons r, and 
temperature T. This represents our equation of state for 
entanglement. 

Vedral [^J, and then Oh and Kim ^(|> nrs f described 
entanglement of noninteracting electron gases. They 
used Green's function approach to find the subsystem of 
the two electrons, separated by a relative distance r. It 
is given by the two-spin density matrix [l(ij |: 

„2 



P12W = -^-[5 ai ,a' l 5 a2 ^ 2 - f(r) Saua'^a'^c 



(1) 



where a± , 02 , a[ , a' 2 are the spin variables of the two elec- 
trons under consideration. Following |S ^| , the state is 
entangled iff the Peres-Horodecki (partial transposition) 
condition tlTj is satisfied : 



f(r) 2 > J , 



(2) 



where / at zero temperature and non zero temperature 
is given by: 

k sinfcr dk (3) 
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and 



f(r,T) 



fkp jq 



k rik sin kr dk , 



(4) 
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where 



(5) 



is the Fermi-Dirac distribution, fi is the chemical poten- 
tial, and £fe are the energy levels. The region of entan- 
glement is: 



< r < r e 



(6) 



where r is the relative distance between the two electrons, 
and the entanglement distance r e is found to be inversely 
proportional to the Fermi momentum kp: 



(7) 



The proportionality constant £ can be found numerically 
by solving f 2 = i. The integral in Eq.QJ can be 
approximated for different temperatures [lOj . In general, 
C is a function of temperature, but for T sufficiently close 
to zero, ( and the entanglement distance r e only change 
slightly from the ones at zero temperature. However, the 
functional dependence of the relation fcj? = — remains 
the same for a gas with constant density. A typical 
value for small temperatures is ( ~ 1.8. We note that 
there will be higher order entanglement present (3, 4 
etc. electrons), but we will not consider it here. 

Eq.@) is the upper bound for the distance of two 
entangled electrons. The entanglement distance r e is of 
the order of the Fermi wavelength, which corresponds 
to the spread of the momentum wavefunction of the 
electrons that are close to the Fermi sphere. In general, 
only these electrons are of physical significance. Eq.© is 
the most restrictive bound for entanglement distance and 
is valid for all electrons in the gas. Any randomly chosen 
electrons that are closer then r e are definitely entangled. 
However, due to the spread in their wavefunction, two 
electrons whose momenta lie far below the Fermi surface 
can also be entangled even if they are further apart 
then r e . This however implies the knowledge of their 
momenta k, which is not accessible in this case. Eq.© 
was derived by making no preselection of momenta. The 
function f(r) in Eq.JQ) is computed by summing over all 
possible momenta, from k = to k — kp, for T = 0. In 
order to have entanglement, f(r ) must then satisfy the 
Peres-Horodecki criterion. In [lflj, the entanglement was 
already discussed as a function of the relative distance 
between the electrons. It was found that the entangle- 
ment is maximal if the relative distance between the 
electrons is zero. In this case, the spatial wavefunctions 
fully overlap, and because of the antisymmetric nature 
of the two electron wavefunction, the particles must 
have opposite spin. No spin configuration (| fj.), | J,|)) is 
discriminated and therefore the two electrons are in the 
maximally entangled Bell state, \*f?~) = Ti) ~ I IT))- 



We now relate this entanglement to P. First we 
discuss the domain of validity of our approach. Let us 
start by considering an electron gas at T = 0. All the 
electrons fill up the quantum states in accordance with 
the Pauli exclusion principle so that the total energy has 
its smallest value. The electrons have momenta from 
p = to pp = hkp = (37r 2 n) 3 h. The radius of the Fermi 
sphere in momentum space is the Fermi momentum 
kF, and n = N/V is the density. The number of states 
in a volume element of width dp is proportional to p 2 . 
The total energy can then be obtained by multiplying 
the number of states by p 2 /2m and integrating over all 
momenta. This leads to the general equation of state 
of a Fermi gas, from which the degeneracy pressure can 
be deduced. We distinguish between the non-relativistic 
and the extreme relativistic case |T^ : 



P 



- — r 1 n 3 non-rel. 

5m 



i — t 713 rel. 



(8) 



where m is the mass of the electron. The non-relativistic 
expression is valid for low temperatures, provided that 
T <C Tp, which is the the Fermi temperature. Rewrit- 



ing this gives T 



2mkf 



-na, which is satisfied if the 



gas density is sufficiently high. A degenerate electron gas 
has the peculiar property that it increasingly approaches 
the ideal gas as the density increases |l2|. This can be 
seen from the following argument. Following Landau and 
Lifshitz |l2| , we consider a plasma of electrons and posi- 
tively charged nuclei. When the electrons move indepen- 
dently as in an ideal gas, the mean kinetic energy (which 
in order of magnitude is equal to the Fermi energy) is 
much higher than the Coulomb interaction. For a nu- 
cleus of charge Ze and with the mean distance between 
the electrons and the nuclei being a: 

Ecoui omb *C Ekinetic 

Ze 2 pp 2 

« ^— 

a 2m 

Ze 2 n^ h 2 n% , ,Z,i 

, — <§; with a~( — ) 3 

2m n 

.e 2 m, , 2 

n » (-jj2-) Z ■ 

There is no electron interaction if the gas density is 
sufficiently high. 

Let us now relate r e to the degeneracy pressure. 
The non-relativistic expression for the degeneracy 
pressure is given by Eq.JHJ, which is valid at low 
temperatures if the density is high. High density means 
no electron interaction, so we use the same treatment 
as in jjj ^J- The extreme relativistic expression for 



3 



the degeneracy pressure is valid at all temperatures. 
The Fermi momentum can be expressed in terms of the 
density n, and Eq.{7J) is rewritten as: 



c 



(3tT 2 ) 



3 77,3 



(9) 



The entanglement distance is approximately equal to 
the average distance between the electrons in the gas. 
Rewriting the degeneracy pressure in terms of r e , we ob- 
tain for the non-rclativistic and extreme rclativistic case: 



P 



15n 2 m T e 



12^ 'e 



" 5 non-rel. 



(10) 



rel. 



This expression is valid if T <C Tp. The constants ( and 
£' are different because in the extreme relativistic limit, 
the density is much higher then in the non-relativistic 
case. When the degeneracy pressure of a given gas is 
high, r e is small and the Fermi momentum is high. This 
means that the electrons are moving at high velocities 
and are localized in space with no overlap of their 
wavefunctions, hence no entanglement is present in this 
limit. 

We now give an estimate for r e in white dwarfs. 
White dwarfs are examples of highly degenerate matter. 
They are objects of extremely high density. We express 
the electron density as: 



Zp 
Am 



H 



(11) 



where Z and A are the number of protons and nucleons 
and rriH is the mass of the Hydrogen atom. If we now 
assume that the density is constant, we can deduce en- 
tanglement distance for a white dwarf of mass M and 
radius R. The density is: 



M 



P : 



4/37ri? 3 
Using Ea. l|ll[l . we write: 

r e oc i?M"3 



(12) 



(13) 



We estimate the entanglement distance of SIRIUS B, 
which is a carbon-oxygen white dwarf, of mass lAf Q , 
radius R = O.OOSi?©, and a temperature of T = 27000fiT 
|l3|. Assuming non-relativistic electrons, and calcu- 
lating the Fermi temperature gives T/Tp ~ 10~ 6 , 
which justifies our low temperature approximation. 
From Eo. H13|l . the entanglement distance is found to 
be r e ~ 6 x lCP 13 m. This is two orders of magnitude 
bigger then the size of a typical nucleus. Note that here 
entanglement length is much smaller than the size R of 



the system. This should be contrasted with the limit of 
entanglement length being of the order of R, i.e r e ~ R, 
for which the mass of the system is m ~ 2.7 x 10~ 27 fcg, 
which is about the mass of a nucleus. 

Let us finally express the entropy of formation as 
a function of the degeneracy pressure, relative distance 
between the electrons r, and temperature T. Relating 
the Fermi momentum to the degeneracy pressure for the 
non-rclativistic and extreme rclativistic case: 



(14) 




rel. 



so Eq.(0J becomes 



rPS 



Jg 00 k rik sin kr dk non-rel. 



Tr,P,T 



(15) 



-^s- f n °° k tit sin kr dk rel. 



with 



ft 2 35 

157r 2 m 



and 7' 



r \ - 



12-k 2 



(16) 



We now describe the system as it is compressed. If 
the density n changes, so does the Fermi momentum 
kp = (37t 2 ti)3 and the degeneracy pressure P. This cor- 
responds to a white dwarf, where the particles are forced 
into a smaller and smaller volume by the gravitational 
pressure. This process can continue until the extreme 
relativistic limit is reached where the energy of the elec- 
trons becomes large compared to mc 2 . 

The integrals in Ea. p5|l can be solved numerically for 
a given temperature. At low temperatures, we take the 
chemical potential to be p, ~ tp. All the information 
needed to compute entanglement is given in the two par- 
ticle density matrix Eq.Q). From this density matrix, 
any entanglement measure E can be computed as a func- 
tion of f r ,p,T- We express the entropy of formation in 
terms of electron distance, degeneracy pressure and tem- 
perature for an electron gas. The equation of state for 
entanglement becomes: 



E F (r,P,T) = h 



1 

5 N 



2 fr,P,T 



(17) 



with h(y) = -y log 2 y - (1 - y) log 2 (l - y), and f r ^ T is 
either the non-relativistic or extreme relativistic expres- 
sion. We can also use the concurrence 14j to express 
entanglement: 



C(r,P,T) 



r,P,T 



1 



fr,P,T 



,0}. 



(18) 
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FIG. 1: Concurrence C (solid line) and Entropy of formation 
Ef (dotted line) at zero temperature with fixed r = 10 _10 m 
as a function of non-relativistic degeneracy pressure P, for a 
gas with initial entanglement distance r e = 10 -8 m. 

Consider a non-relativistic gas with high density at zero 
temperature, which undergoes compression. This is 
shown in Fig. 1. If we fix the relative distance r of the 
two electrons to be smaller then the entanglement dis- 
tance r e , then the two electrons are certainly entangled. 
The smaller the distance r is, the more entangled are the 
electrons. If we now compress the system, the density n 
increases, and r e in Eq.@ decreases. The entanglement 
distance r e approaches the relative electron distance 
r from above, and the entanglement between the two 
electrons decreases. If the system is compressed further, 
r e becomes smaller then r, and entanglement vanishes. 
This means that the electrons are forced into a smaller 
and smaller volume, so their momenta increase and their 
wavefunctions overlap less, until their is no overlap and 
therefore no entanglement. 

Ea. (|17|l or Eq. (|18|l are not equations of state in 
the " classical sense" of that word. We call them like that 
because they treat entanglement on the same footing 
as pressure and temperature, which are classical ther- 
modynamical variables. The only microscopic variable 
left is the relative distance between the electrons. If 
we do not want a microscopic variable to appear in the 
equation of state, we can average over the region where 
entanglement between two electrons certainly exists. 
For convenience, first define the variable x = UfT. The 
average then implies integrating from x = to x = x{T). 
We find x(T) by solving f(x,T) 2 — |. We can now 
find the average amount of entanglement for electron 
distances where there exists entanglement, using a 



general integrable entanglement measure E: 

, r x=x(T) 

{E{T)) = —- E(x,T)dx. (19) 

xy 1 ) Jx=o 



In conclusion, we have considered a degenerate electron 
gas with high density at low temperatures, and an 
extremely relativistic electron gas. We first related 
entanglement distance to degeneracy pressure. We then 
estimated entanglement distance for the simplest model 
of a white dwarf. For the first time, the entropy of for- 
mation and the concurrence was related to temperature, 
pressure and relative distance between the electrons. 
It is valid for temperatures which are sufficiently close 
to zero and high density for the non-relativistic case, 
and for all temperatures for the extreme relativistic 
case. In this work entanglement is purely due to particle 
statistics (Pauli exclusion principle) and not to fermion 
interaction. Future work would consist of investigating 
entanglement for fermions that physically interact. It 
would also be interesting to extend the analysis to 
multipartite entanglement and to see if one could indeed 
treat entanglement as a thermodynamical variable, and 
write an equation for entanglement similar to the the 
fundamental thermodynamic equation for energy U, 
dU = TdS - PdV. 

We are grateful to Marcelo Franga Santos for helpful 
discussions. 



* Electronic address: christian.lunkes@ic.ac.uk 
' Electronic address: c.brukner@ic.ac.uk 

* Electronic address: v.vedral@ic.ac.uk 

[1] V. Vedral and E. Kashefi, Phys. Rev. Lett. 89, 037903, 
(2002). 

[2] M. B. Plenio and V. Vedral, Contemp. Phys. 39, 431, 
(1998). 

[3] V. Vedral, New J. Phys. 6, 22, (2004). 
[4] S. Popescu and D. Rohrlich, Phys. Rev. A 56, R3319 
(1997). 

[5] M. C. Arnesen, S. Bose, V. Vedral, Phys. Rev. Lett. 87, 
017901 (2001). 

[6] S. Gosh, T. F. Rosenbaum, G. Aeppli, S. N. Copper- 
smith, Nature 425, 48 (2003). 

[7] C. H. Bennett et al, Phys. Rev. A 54, 3824 (1996). 

[8] C. Brukner and V. Vedral, e-print |quant-p h/0406040 
(2004). 

[9] V. Vedral, CEJP 2, 289 (2003). 
[10] S. Oh and J. Kim, Phys. Rev. A 69, 054305 (2004). 
[11] A. Peres, Phys. Lett. A 202, 16 (1995). 
[12] Landau and Lifshitz, Statistical Physics 1, (Butterworth- 

Heinemann, 1980). 
[13] B. W. Carroll and D. A. Ostlie, Modern Astrophysics, 

(Addison- Wesley, 1996) . 
[14] V. Coffman, J. Kundu, W. Wootters, Phys. Rev. A 61, 

052306 (2000). 



